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I.  INTRODUCTION 

When  diffuse  objects  are  illuminated  with  coherent 
radiation,  the  appearance  of  the  image  is  degraded  by  the 
presence  of  speckle.  The  effect  of  the  speckle  is  to  super¬ 
impose  a  noise-like  structure  which  masks  the  spatial  infor¬ 
mation  present  in  the  image. 

Considerable  work  has  been  done  in  characterizing  the 
speckle  over  the  last  few  years  [1] .  Both  theoretical  and 
experimental  work  has  revealed  new  ideas  which  have  enhanced 
our  understanding  of  the  speckle  phenomena. 

The  purpose  of  this  work  is  to  extend  the  statistical 
techniques  of  hypothesis  testing  first  used  by  Dainty  [2]  to 
study  speckle  and  to  apply  these  extended  techniques  to  the 
data  generated  in  the  studies  conducted  by  Guenther  et  al. 
[3,4]. 


Dainty  used  simple  statistical  detection  theory  to 
study  the  problem  of  detecting  small  images  immersed  in  a 
background  of  laser-produced  speckle.  He  assumed  that  the 
detection  device  was  a  flying-spot  scanner  and  showed  that 
the  size  and  transmittance  of  the  scanning  aperture  was 
best  determined  by  statistical  considerations  and  not  by 
signal-to-noise  (SNR)  criteria.  He  also  showed  that  the 
typical  detection  problem  of  locating  opaque  disks  in  a 
speckle  background  requires  that  the  geometrical  image  diam¬ 
eter  be  at  least  four  times  the  Rayleigh  resolution  limit 
of  the  lens  used  to  form  the  image  for  reliable  detection. 
The  work  by  Dainty  was  theoretical;  no  experimental  data 
were  presented  to  test  his  theory. 

On  the  other  hand,  the  work  in  References  3  and  4 
was  largely  experimental  with  estimates  of  detection  deter¬ 
mined  empirically.  In  that  work,  a  special  test  target 
patterned  after  one  first  used  by  Rose  is  used.  The  target 
has  disks  of  varying  transmissivity  and  diameter  which 
are  imaged  in  the  presence  of  speckle.  The  transmittance 
of  the  disks  varies  from  0.88  to  0.012,  while  the  size 
varies  from  about  two  times  the  speckle  size  to  about  14 
times  the  speckle  size.  The  test  target  is  imaged  under 
various  illumination  conditions. 

In  the  following  we  will  present,  in  a  clear  way,  the 
theory  first  proposed  by  Dainty  and  show  how  this  theory  can 
apply  to  the  experimental  work  of  Guenther  et  al. 
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II.  TESTING  OF  SIMPLE  HYPOTHESIS  [5] 

In  hypothesis  testing,  an  observer  measures  a  single 
quantity  X  and  on  the  basis  of  this  measurement  chooses 
between  two  hypotheses  Hj  or  Hp  If  X  always  was  equal  to 
ag  when  the  hypothesis  Hg  was  true  and  to'  a^  when  was 
true,  there  would  be  no  need  for  statistical  hypothesis 
testing.  However,  because  we  may  make  errors  in  measure¬ 
ment  or  we  have  factors  operating  which  we  do  not  under¬ 
stand,  the  observation  quantity  X  is  a  random  variable  that 
we  can  describe  statistically  by  giving  a  probability 
density  function  pg (X)  or  p^(X)  if  Hg  or  Hi  applies.  The 
hypotheses  are  called  simple  if  both  pg(X)  and  pi (X)  are 
completely  known  and  they  do  not  depend  on  unknown  para¬ 
meters  . 


A  strategy  must  be  adopted  by  the  observer  which 
assigns  one  hypothesis  or  the  other  to  each  observable  out¬ 
come.  This  strategy  divides  the  possible  values  of  X  into 
two  ranges,  Rg  and  R^ ,  so  that  if  X  lies  in  Rg ,  Hg  is 
chosen  and  if  X  lies  in  Rj_,  H^  is  chosen.  The  regions 
should  be  so  determined  as  to  get  the  best  average  results 
in  the  experiment. 

A  simple  example  of  this  statistical  test  is  to 
choose  between  two  hypotheses: 


Ho:  rain  tomorrow 

Hy:  no  rain  tomorrow 

A  possible  way  of  obtaining  such  a  decision  is  to  determine 
the  average  rate  of  change  of  barometric  pressure  during 
the  past  24  hours. 

Suppose  that 

Pg (X)  is  the  p.d.f.  of  rates  X  on  days  before  rain  = 

2  2l 

ag)  /2a 

of  rates  X  on  days  before  no 


2  ira 


— 


exp 


(X  - 


Pi  (X)  is  the  p.d.f. 


rain  = 


2tto" 


-U 


exp  jj  (X  -  a^)  2/2o2J 


where  ag  <  0  <  ai , 
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Figure  1.  Probability  density  functions  under  hypotheses 
Hq  and  Hp. 


These  p.d.f.'s  are  shown  in  Figure  1 .  We  can  see  that  on 
days  before  rain  the  barometer  fell  with  an  average  rate 
an  which  is  negative,  while  on  days  before  no  rain  the  aver 
age  rate  is  positive  because  the  barometer  was  rising.  In 
addition  to  the  above  data  contained  in  the  p.d.f. ,  suppose 
we  know  the  fraction  t;  of  rainy  days  that  occurred  in  the 
past.  This  number  x,  can  be  called  the  prior  probability  of 
rain  and  (1  -  5)  the  prior  probability  of  no  rain. 

As  we  can  see  in  Figure  1,  the  p.d.f. 's  overlap; 
thus,  no  matter  how  we  choose  regions  Rq  and  R]_,  we  will 
occasionally  make  a  wrong  decision.  In  fact,  the  probabil 
ity  of  choosing  Hi  when  Ho  is  true,  so-called  error  of  the 
first  kind,  is 


Qq  =  /  P0(x>  dx' 


shown  as  the  crosshatched  area  under  the  po(X)  curve.  The 
probability  of  choosing  Hq  when  Hi  is  true,  an  error  of  the 
second  kind,  is 
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The  value  of  Xo  which  the  observer  will  choose  depends  on 
how  much  wrong  decisions  cost.  Support  that  Cq  is  the  cost 
of  an  error  of  the  first  kind,  Ci  is  the  cost  of  an  error  of 
the  second  kind,  and  CqQo  and  C1Q1  are  called  the  risk  asso¬ 
ciated  with  Hg  and  H]_,  respectively.  Then  the  average  risk 
of  each  decision  is 


C(X0)  =  cCgQg  +  (1  -  e)  CiQi 

CO  X0 

=  cCg  /  pg  (X)  dx  +  ( 1  —  ^  )  C]_  /  Pi  (X)  dx. 

XQ 

We  can  minimize  the  above  by  differentiating  with  respect  to 
Xg  and  putting  the  result  equal  to  zero.  Thus, 


Pi  (Xp)  _  <C0 _ 

PgTXg)  (1  -~OCX  ~  A0 
which  yields 

Xg  =  (ag  +  a\)  /2  +  ai  --o  m  (1-!Cj>-)  Cl  • 

Note  that  if  ^  =  1/2  and  Cg  =  C^,  we  get  Xg  =  (ag  +  a^)/2. 


III.  BAYES  CRITERION 

The  Bayes  criterion  is  a  rule  which  allows  the  deci¬ 
sion  strategy  to  be  picked  so  as  to  minimize  the  average 
risk. 

There  are  four  types  of  costs  possible  which  can  be 
arrayed  in  a  cost  matrix  C. 


where  CjH  is  the  cost  of  choosing  when  Hj  is  true  (i,  j  = 
0,  1).  The  relative  cost  of  error  of  the  first  kind  is  C^g  - 
Cgg  and  of  the  second  kind  Cqi  -  Ci]_,  with  both  relative 
costs  being  positive.  In  the  example. 


C  = 


In  the  general  case  where  elements  of  the  cost  matrix 
are  not  zero,  we  have  the  average  per  decision  as 


C  =  5 


0 


'00 


/  p  (X)  dx  +  C  n  /  p  (X)  dx ; 


+  (1  -  C> 


X 


10  Y 

xo 


0 


C01  /  Pl(X)  dx  +  C11  /  Pj(X)  dx  . 
-  —  X0 


It  can  be  shown  that  the  average  risk  is  minimum  when 
Xg  is  chosen  so  that,  given  the  likelihood  ratio 


A  (X)  =  Pi  (X)  /p0  (X)  , 

Rq  consists  of  values  of  X  for  which  A(X)  <  Ag  and  con¬ 
sists  of  values  of  X  for  which  A(X)  >  Aq  where  Aq  is 


C(C 


Ao  = 


10 


coo) 


(l  -  O  (c 


01 


-  cll} 


If  5,  the  prior  probability  is  known  and  values  can  be 
assigned  to  Cji  and  the  threshold  Xg  dividing  regions  Rg 
and  can  be  found  which  minimizes  the  average  risk  for 
each  decision.  For  the  general  cost  of  the  detection  of 
signals  in  noise,  it  is  often  the  case  that  we  cannot 
reasonably  assign  risks  nor  can  we  know  prior  probabilities; 
the  position  of  the  threshold  is  set  from  experiment  by 
trial  and  error. 


IV.  SIGNALS  IN  NOISE 

The  theory  outlined  above  is  often  used  in  the 
design  of  communications  and  radar  systems.  There,  often 
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the  problem  is  to  determine  whether  a  signal  and  additive 
noise  are  being  observed  or  whether  noise  alone  is  present. 


In  this  case,  the  situation  is  illustrated  in  Figure 
2.  Here  the  crosshatched  area  under  the  p.d.f.  curve  of 
signal  +  noise  gives  the  probability  of  detection  while  the 
doubly  crosshatched  area  A2  under  the  p.d.f.  curve  of  noise 
alone  is  the  probability  of  a  false  alarm.  The  problem  here 
is  determining  the  proper  signal- to-noise  ratio  and  threshold 
to  give  an  adequate  probability  of  detection  while  keeping 
the  probability  of  a  false  alarm  properly  low.  The  SNR, 
which  determines  the  relative  position  of  the  p.d.f.  of 
signal  +  noise  to  the  p.d.f.  of  noise  alone,  along  with  the 
proper  choice  of  threshold,  allows  obtaining  proper  values 
of  detection  and  false  alarms. 


A^=probability  of  detection  A2=probability  of  false  alarm 

Figure  2.  Probability  of  noise  alone  and  signal  +  noise 

illustrating  the  process  of  threshold  detection. 


V.  DISK  IMAGES  IN  SPECKLE 

Suppose  we  have  disk  objects  of  intensity  trans¬ 
missivity  <Ti>  backed  by  a  ground  glass  with  average  trans¬ 
missivity  of  <T_ >  where  we  assume  <Tq>  >  <T.>.  The  diam¬ 
eter  of  the  disk  is  given  by  d  wherea  >  d^  The  quantity 


ds  is  the  approximate  size  of  the  speckle  as  determined  by 
the  Rayleigh  criterion.  That  is,  dg  =  -yy,  where  F  is  the 
focal  length  of  the  lens  and  X  is  the  wavelength  of  the 
light.  Figure  3  illustrates  this  situation. 

Suppose  further  that  we  place  a  photocell  at  plane  P 
with  a  diameter  of  dp.  This  photocell  is  used  to  measure 
the  amount  of  light  which  is  present  at  this  plane.  We  can 
use  the  voltage  reading  from  the  photocell  to  determine  if 
a  disk  is  present  or  not.  The  intensity  measured  by  the 
photocell  is 


where 


I  =  //  I (x ,y)  B(x,y)  dxdy 

^  —  oo 

I(x,y)  is  the  intensity  distribution  at  P3 


and 


(1) 


B(x,y)  =  1  for  x2  +  y2  £  dp  (2) 
=  0  otherwise; 

I(x,y)  can  either  be  the  intensity  due  to  the  presence  of 
a  disk  in  which  case  its  average  value  as  a  function  of  x 
and  y  is  given  by 


<1 (x,y) >d  =  |0A(x,y)|2  *  | K (x , y ) | 2 


(3) 


where 

I  0A  (x , y )  |  2  =  I  <Ti>  < Tg>  for  x2  +  y2  d2  (4) 

=  I< Tg>  otherwise; 

I  is  the  uniform  intensity  of  the  laser  illumination  and 
|K(x,y)|2  is  the  intensity  point  spread  function  of  the 
optical  system  given  by 


K(x,y) 


kD‘ 


8F 


2 

r  j 

^kDrn 

0  1 

L  Ji  1 

l  2F  1 

kDrn 

1 _  _ 

-J 

2F 

(5) 


where  k  =  2ir/X,  J^(X)  is  the  first-order  Bessel  function 
and  r$  =  x^  +  y2.  I(x,y)  will  vary  around  <I(x,y)>  due  to 
the  presence  of  speckle  superimposed  over  the  disk  image  by 
the  ground  glass. 
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GROUND  GLASS 


Figure  3.  Imaging  a  dark  disk  with  a  coherent  optical  sys¬ 
tem  in  the  presence  of  speckle  noise. 


If  there  is  no  disk  present  in  plane  P]_,  then  the 
intensity  I(x,y)  is  the  intensity  of  the  imaged  ground  glass. 
The  average  value  of  this  intensity  is  given  by 

<Ilx,y)>nd  =  I<V-  (6) 

Again,  I(x,y)  will  vary  around  the  average  -I-  because  of 
the  ground  glass. 

The  p.d.f.  associated  with  the  photocell  reading  due 
to  the  intensity  variation  when  no  disk  is  present  is  given 
by  Pq(I),  while  the  p.d.f.  associated  with  the  photocell 
reading  when  a  disk  is  present  is  given  by  p g ( 1) .  The 
p.d.f.'s  are  shown  in  Figure  4.  Generally,  the  p.d.f.'s 
will  overlap  and  with  a  threshold  I  sometimes  the  wrong 
decision  will  be  made  whether  a  disk  is  or  is  not  present. 

Two  kinds  of  errors  can  be  made.  Errors  of  the  first  kind 
are  the  apparent  detection  of  an  image  when  it  is  not  pre¬ 
sent.  This  false  alarm  probably  a  is  given  by 
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and  the  distance  between  the  p.d.f.'s  is  determined  by  the 
average  value  of  the  intensity  transmittance  of  the  disk 
and  ground  glass. 


VI.  PROBABILITY  DENSITY  FUNCTIONS 


An  approximate  expression  for  Pq  and  Pi  [6]  is 


for  I 
tion , 
is 


>  0,  zero  otherwise.  Where 
M  =  (<I>/oT) 2 ,  <I>  average 
the  standard  deviation. 


r  (x) 

value 


is 

of 


the  gamma  func- 
the  p . d . f .  and 


VII.  CALCULATION  OF  a  AND  (3 

In  order  to  proceed  with  the  calculation  of  •  and  ,•  , 
we  need  to  calculate  the  mean  <I>  and  the  standard  devia¬ 
tion  oj  of  the  p.d.f.'s  when  an  image  is  present  and  not 
present. 

pn(I)  [No  Disk  Present] 


(0X>nd  =  <J>nd  H  Nnd(u'  v)  'b(u'  v)  I  2  dudv 

—  oo 

where  b(u,  v)  is  the  Fourier  transform  of  B(x,  y) . 

Nnd(u'  v)  =  KndlF(u'  v)!2  *  If(Tu'  “v)i2 

where  F(u,  v)  is  the  aperture  of  the  filter  in  plane  P2 
The  constant  Knd  can  be  evaluated  from  the  fact  that 


(10) 

(ID 


oo 

/  Po  ( I )  dl  =  1.  (13) 

—  oo 

Once  <I>nd  and  (°i)nd  are  evaluated,  they  can  be  substituted 
into  Equation  (9)  to  obtain  p0(I). 


12 


Pi ( I )  [Disk  Present] 


<I>d  is  obtained  through  the  evaluation  of  Equation 
(3)  using  Equations  (4)  and  (5),  and  (oj)d  is  given  by 

(Oj)^  =  <I>d  //  Nd(u,  v)  ! b ( u ,  v) !  dudv  (14) 

—  oo 

where  b(u,  v)  is  again  the  Fourier  transform  of  B(x,  y)  and 

N , (u ,  v)  =  K,|F(u,  v)  I2  *  j  F (-u ,  -v)  !2  (15) 

d  d 

where  F  is  the  same  as  above.  Kd  is  evaluated  as  above  by 
noting  that 
00 

/  P-^I)  dl  =  1.  (16) 

—  oo 

When  the  two  p.d.f.'s  are  found,  a  and  2  can  be  calculated 


a 


P0(I)  dl 


a 


M 


<I> 


nd 


M 


tM-1 

I  exp 


r  (M) 


MI 


<I> 


nd 


dl 


where 


M 


<I> 


(a. 


2 

nd 

2 

nd 


IT 

6  =  /  p,  (I)  dl 
0 


=  / 
0 


M 


<I> 


M 


tM-1 

I  exp 


r  (M) 


MI 


<I> 


dl 


(17) 


(18) 
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where 


<I> 


M  = 


<U>a 


VIII.  CONCLUSIONS 


We  can  note 

how  the 

size  of  the  disk, 

the 

intensity 

transmissivity  and 

the  size 

of  the  scanning 

aperture  affect 

the  calculation  of 

i  and 

Parameter 

Ol 

- 

d 

size  of  disk 

Eq . 

(3)  *!,d 

Ti 

disk  transmittance 

- 

Eq. 

(3)  -I  d 

Ts 

g.g.  transmittance 

Eq . 

(6)  'r>nd 

Eq. 

(3)  Id 

D 

size  of  aperture 

Eq . 

<6>  *  nd 

Eq. 

(3)  Id 

in  plane  P2 

Eq . 

lil>  <°I>nd 

Eq . 

(14)  (ox)d 

dP 

size  of  scanning 

Eq. 

<X>nd 

Eq . 

(1)  I>d 

aperture 

Eq. 

(11>  <Vnd 

Eq . 

(14)  (cI)d 
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